Exact non-equilibrium current from the partition function 
for impurity transport problems 
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We study the partition functions of quantum impurity problems in the domain of complex a pplied 



bias for its re lation to the non-equilibrium current suggested by Fendley, Lesage and Saleur (:ond- 
mat/9510055). The problem is reformulated as a certain generalization of the linear response theory 
that accomodates an additional complex variable. It is shown that the mentioned relation holds 
in a rather generic case in the linear response limit, or under certain condition out of equilibrium. 
This condition is trivially satisfied by the quadratic Hamiltonians and is rather restrictive for the 
interacting models. An example is given when the condition is violated. 

PACS numbers; 72.15.Qm, 73.40. Gk, 72.10.Bg, 71.30.+h 



I. INTRODUCTION 

Strongly correlated quantum systems out of equilib- 
rium is one of the lesser understood domains in theoret- 
ical physics, mainly because such systems reveal much 
more complex behavior compared to the case of equi- 
librium. Correspondinly, there are not that many non- 
trivial examples known today that are solved exactly out 
of equilibrium, and extension of powerful field-theoretical 
methods to non-equilibrium situations is desirable. 

Situation is somewhat simpler in the stationary non- 
equilibrium regime. In a remarkable work, Ref. the 
authors put forward a novel idea of how to calculate 
the current for integrable systems out of equilibrium. In 
particular, they obtained, based on the standard Bethe 
ansatz technique, a non-perturbative expression for the 
backscattering current through a single impurity in the 
quantum Hall bar with v ^ ^. Their result, cast into 
appropriate form H, reads: 
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where T stands for temperature, = ^ is a filling frac- 
tion, Xbs is an effective tunneling strength of the im- 
purity, 2imp is impurity's partition function defined as 
the ratio of the equilibrium partition function of the full 
system to that of the "pure" system without impurity 
(the leads), and analytic continuation from the discrete 
imaginary to real chemical potential of the edges fi is as- 
sumed. It is worth emphasizing that Eq. (Q) was meant 
to have a truly non-equilibrium nature, unlike similar ex- 
pressions for persistent currents |^,^, where the current 
/ = cdF/d^ is a purely equilibrium phenomenon, a con- 
sequence of the incident flux $ on the structure of energy 
levels. 

Relation (^ was obtained through several intermedi- 
ate steps, rather than derived directly, (a) First, the "fu- 
sion" of Bethe ansatz with the Bolzmann rate equation 



was implemented as follows The quantum Hall bar 
with impurity, described by the Luttinger theory with 
point-hke backscattering term XbsS{x){'if'l'^ r + ^f^^fi) 
was reformulated in terms of the boundary sine-Gordon 
model 
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which is exactly solvable in equilibrium and can be diago- 
nalized by the Bethe ansatz. Transport of Laughlin (e/3) 
quasiparticles through impurity Q maps onto a scatter- 
ing of sine-Gordon quasiparticles (kinks, antikinks and 
breathers) off the boundary. All the interaction in (||) 
is now at the boundary, which behaves like a non-elastic 
scatterer: kink can be bounced as an antikink and vice 
versa, which changes the charge of the system. However, 
it is a bare Hamiltonian where the bulk interactions are 
absent. As a result of diagonalization, sine-Gordon quasi- 
particles interact also in the bulk by a pairwise point-like 
interaction that adds merely a phase-shift with momen- 
tum of each quasiparticle preserved - a consequence of 
the peculiar conservation laws of (|^), and the distribution 
function of the gas of quasiparticles differs from the usual 
Fermi one. Under these circumstances, the current was 
calculated on the basis of the standard probability argu- 
ments which are known to lead to the (classical) Bolz- 
mann rate equation [Q : 

e f°° 

Ib = J dp[n+{p,fi,T)-n_{p,fi,T)]\S^{p,Xbs)\'^- 

(3) 

This expression was referred in literature to as an ex- 
act result because of the exact quantum nature of all the 
entries of expression (^ used in Q. In particular, the 
densities of states n± of quasiparticles in the presence 
of the external bias /i dx^ were found from the ther- 
modynamic Bethe ansatz (TBA) 0,0], whereas for the 
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transition probability W = |(+|— )P out of equilibrium 
it was taken the kink-antikink equilibrium scattering ma- 
trix element , obtained in [|j . 

(b) Second, the exact impurity partition function of (|^) 
was obtained either as a Coulomb gas expansion, or based 
on the relation to the Kondo model in a magnetic field 
and using TBA for the latter |||. Eventually, formula (0) 
was checked by comparing series expansions of both sides 
to a high order HJ^. 

Existence of valid generic relations like Eq. (|l|) would 
be rather helpful. Indeed, a partition function is much 
easier to find and can be tackled by a variety of tech- 
niques: Feynmann integration in imaginary time, Monte 
Carlo simulations, Bethe ansatz etc. Then the proper an- 
alytic continuation in /Lt must be done. Speaking loosely, 
it is in some sense analogous to the analytic continua- 
tion from Matsubara frequencies to the real axis in the 
standard linear response theory. Out of equilibrium, the 
real time retarded Green function has no obvious thermo- 
dynamic interpretation on the imaginary frequency axis; 
equations of the kind (||) offer such an interpretation. 
Eq. allows to cook up the non-equilibrium retarded 
Green function in the complex frequency plane from the 
Matrsubara Green functions in a very simple way. The 
important question is, however, how general formula (|^) 
is. 

In this paper we undertake further studies of the re- 
lations of the form of Eq. (|l|) without appealing to the 
knowledge of Bcthe ansatz, but rather on the basis of 
standard Green functions technique. We focuse on the 
transport properties of one-dimensional electrons pass- 
ing through an impurity or an interacting region, the 
quantity of interest being an I-V characteristic. Exper- 
imentally, such systems are realized, e. g. , as quantum 
dots, quenched nanoconstrictions, quantum wires with 
an impurity or the edges of quantum Hall liquid. On 
the theoretical level, we start with a hypothetical model 
of two reservoirs, L and R, of free ID electrons with a 
linearized spectrum, kept at different constant chemical 
potentials fi^ a-nd /i/j, and described by the Hamiltonian 
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(4) 



where we have chosen left-moving branch for both L and 
R species and set = 1 (similar starting point 
of view was taken in Q|, the reservoirs being the edge 
channels or kinks-antikinks in the sine-Gordon formula- 
tion). Each of the reservoirs (referred to as leads below) 
is intially in a well-defined equilibrium state. Then, one 
allows for a current to flow from L to i? by adiabati- 
cally switching on some coupling term. Since the amount 
of electrons in reservoirs is infinite in the infinite band- 
width limit, the reservoirs never equilibrate and the sys- 
tem achieves some non-equilibrium steady state with a 
constant flux in the energy space. Equation for an exact 



quantum current can be derived rigorously in the spirit 

offnl. 

To be specific, we consider three types of coupling 
terms: point-contact, resonant-level and Kondo. Inter- 
actions can be introduced in the free (quadratic) mod- 
els (a) by moving away from the Toulouse point in the 
Kondo model; (b) by taking leads to be the Luttinger liq- 
uid instead of free electrons in the point-contact model 
and (c) by considering a larger space of interacting states 
{dn} instead of one state in the resonant-level model. 
We study in detail only the latter case here. We argue 
that: a) equation (|l|) holds for free systems (quadratic 
couphngs), and b) Eq. (|l|) holds for the interacting sys- 
tems under certain condition for the self-energy, or as a 
generic formula in the linear response limit. The results 
can be easily generilized to the case of transport through 
an interacting region if one interprets {dn\ as a set of 
eigenstates in the isolated region labeled by the quan- 
tum number n with the hopping strength depending on 



II. PARTITION FUNCTION vs. 
NON-EQUILIBRIUM CURRENT 

Consider the following Hamiltonian, which describes 
the leads ('0), the impurity {d) located at x = 0, and 
possible interactions {Hint)'- 

H = Ho + \Y, (^/'+(0)rf + d+^m(0)) 

771 — L.R 

+ eod+d + H,„t{{d}), (5) 

Operator c?"*" creates an electron on the impurity site at 
the energy cq. For this model, an analogue of relation (|l|) 
to be shown reads: 



t7np 



where F = Q. 

Define the impurity partition function as 



(6) 



rY^.^-(H+itj.LNL+itJ.RNR)/T 



(7) 



where ifJ,L{R) take the values on the bosonic Matsubara 
frequencies. The key point in Eq. (^) is analytic contin- 
uation in the chemical potential. To perform it carefully, 
we proceed as follows. From the definition of Zimp and 
Eq. (0) one has: 
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— log Z,mp = -T-^ (V^™'^ + 
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where G^^^ = - j^''^ dre^'^-^ {Trd{T)d+ (Q)) is the full 
Matsubara propagator for the impurity, 



ydd 



eo - TiGoL + Gor) - J:^{iujn\{ifi}) ' 



(10) 



LtJn are the fermionic Matsubara frequencies, and to ob- 
tain (Q) the Dyson equations were used, which follow 
from the quadratic nature of coupling between the impu- 
rity and the leads: 



{Tr^Ld+),=o = -XGoL G^d, 

{Trdlpl)a:=0 = ->^GQRGdd, 
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(lib) 



where 
GoLiR) 
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are free Matrubara Green functions of the disconnected 
leads (A = 0) at X = 0. Substituting |l3] 



GoL{R)i^n) = lim / ^ 
A^oo 27r 
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_A ^TT lUJn+p- ^^^'L(R) 

-^sign(w„ - AiL(i?,)), 



(13) 



into (M) one has: 



— logZ„„p = i ^ [GddiiuJn + i max[^L, /^i?]) 
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Qdd{-^^n + iTm■n[^iL,^J.R])\ , (14) 



where it was used that ifJ.L{R) ^re restricted to the bosonic 
Matsubara frequencies, 2TriNT, which allows to bring it 
in the argument of the Green functions. Emergence of the 
non-analytic functions "max" and "min" at this stage of 
calculation is an important feature. Expression ( p^ can 
be viewed as a function of two complex variables iil and 
/ifl, with a cut at — Mfl llll- Because of the cut, one 
can define two separate analytic continuations: once from 
the domain > jj,R, and once from > to get two 
different functions, related to and Z~^p (this pro- 

cedure is analogous to obtaining retarded and advanced 
Green functions from the Matsubara Green function by 
continuation in lu). Subtracting one function from an- 
other and substituting ii-J-L(R) ~* P'-l{r)i one obtains: 



d Z^ 

log = i ^ [Qddii^n + ml) - Gddii^n + Mi?) 

imp ij„>0 

+ Gddi-i^n + ml) - Gddi-i^n + Mfl)] , (15) 



dT 



where the continuation ^ fi is implicitly assumed in 
( p^ ) also in the second argument, {/i}, which is omitted 
here for brievity, and the expression in the square brack- 
ets contains now only analytic functions over real fi. Note 



that the sum over w„ in (|j) diverges and hence an ana- 
lytic continuation in /i might be not unique. We assume 
the sum in (^ to be regularized in the standard way by 
inserting exp(ia;„r) ll^. Then, since expression in the 
right-hand side of (|14[) is now a function of /i which van- 
ishes when oo (by the reason similar to as a Green 
function vanishes when iujn oo), the analytic contin- 
uation in fi to the appropriate domain in the complex fi 
plane becomes unique by the well-known mathematical 
theorem. One may also notice that, without the prelimi- 
nary regularization, the divergence in Z'^ cancels anainst 
the same divergence in Z~ and the sum in (|l5|) is conver- 
gent, which also defines a certain way of regularization. 
The correctness of this procedure can be checked against 
the non-interacting case. 

Next, consider an interacting non-equilibrium system 
in the real time. By definition. 



J dpL i ITT 1 



(16) 



where pl = sfpi^fpL is the density of charge in the left 
lead. Opening the commutator in ( |l6| ) and taking the 
average, one obtains 



Jl 



ieX 



{{d+{t)Mt,0))-{^tit,0)dm. (17) 



Using Dyson equations (Eq. (B13) from appendix B), in 



analogy with the Matsubara calculation above, correla- 
tors in (17) can be expressed in terms of the free lead 
propagators and the full c?-level retarded Green function 
to obtain ||ll|,|l|l: 



J = h{JL — Jr) 



~h 



dio[fL~fR]lraGfd{u^), (18) 



where /l.r = /(<^ — Hl,r) are the Fermi functions. Eq. 
( |T8| ) was derived without any reference to the nature of 
impurity interactions Hint] in the case when d has ad- 
ditional quantum numbers, trace operation must be per- 
formed |l|. Full d — d propagator. 



riR — 

'^dd — 
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can be obtained by the Keldysh technique if one uses 
( |l9| ) with = as a zeroth-order approximation for the 
impurity Green function and then adiabatically switches 
on interactions. Using that 2iImG^ — — G^, and 
that G^ (G'^) is analytic in the upper (lower) complex 
half-plane, we split the expression for current ( [T^ ) into 
two parts, one containing G^ and another G'^ and close 
the contour in the upper/lower complex half-planes re- 
spectively. Then, by the residue theorem, one has: 
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+ Gdd{^^B^ - i^n) - G^^ipL - icj„)] 



(20) 



3 



To this formula is related our first observation: find- 
ing the exact quantum current is equivalent to evalu- 
ating the non-equlibrium retarded Green function in the 
complex Lo plane at the fermionic Matsubara frequencies 
shifted from the imaginary axis by the amount of the 
chemical potential. Now it is convenient to compare the 
real time calculation with the Matsubara claculation, Eq. 
( p^ ) , since both expressions are defined over the same set 
of points in the complex lo plane {fJ-L{R) + i^n)- For the 
non-interacting case self-energy vanishes, so that one has 
an exact expression 



1 



UJ — Eq + ir 



(21) 



Note that in the absense of interactions the retarded non- 
equilibrium Green function does not depend on the chem- 
ical potentials of the leads and coincides with the equi- 
librium retarded Green function. With interactions, the 
self-energy (and thus G^j) appears to be in general bias- 
dependent. For the current one also has an exact result, 
given in Appendix A. Comparing ( p^ with ( |20| ) for the 
free correlators (S*^ — = 0) we see that they agree 
term by term, hence Eq. (^ holds. We checked also that 
the analogue of Eq. (^) holds for other free models, e.g. 
the anisotropic two-channel Kondo model in the Toulouse 
limit @jl9| or the point-contact model (see Appendix 
A): 



Hpc — Hq 



(22) 



Hamiltonian ( |22| ) corresponds to a particular case of the 
quantum Hall bar with an impurity ||^ a,t v = 1, when 
the Luttinger liquid on edges reduces to free fermions, 
and thus it provides a direct check of Eq. (|l|). However, 
in order to satisfy Eq. (^) in the presence of interactions, 
one needs, e. g., 

S^(zc^„ + - S*^(ic^„ + iAiL|{«/i})v-M (23) 

for ujn > 0, or, equivalently. 



(24) 



Equation (^3|) states that, out of equilibrium, if the re- 
tarded self-energy is equal to the Matsubara self-energy 
analytically continued both in uj and in /i, then Eq. (|) 
holds, and it is very likely that for generic interactions 
this statement works both ways. Thus, Eq. (^3|) repre- 
sents a very interesting conjecture. It is a very simple 
generalization of the analagous result in the linear re- 
sponse theory. Indeed, to obtain the current in the linear 
response limit, //L,/i_R — > [Q, one can substitute for 
the retarded Green function in Eq. (p^ ) its equilibrium 
value, and then Eq. ( p3|) becomes just a familiar relation 
from the linear response theory. Hence, Eq. (||) holds in 
the linear response limit. 



III. PERTURBATIVE ANALYSIS OUT OF 
EQUILIBRIUM 

In strongly interacting ID systems a naive perturba- 
tion series without some self-consistent resummations of- 
ten do not reveal the correct physical properties of the 
system. However, such a series, if convergent, are quite 
suitable for the purposes of comparison of the mathemat- 
ical structure of both sides of ( p3| ) . 

To see that in general out of equilibrium Eq. ( |23| ) is not 
satisfied, consider the simplest interacting model with 
{dn} = {d-f ) '^il ^-iid interaction of the form 



Hint = Ud^I'd^'d^di, 



(25) 



(correspondingly, electrons in the leads also carry spin 
and interactions are such that the total spin is preserved). 
This is known as the Anderson model pl[ , and it was 
studied in the U = oo p^-[24|] as well as small U/T limits 
p5[ out of equilibrium for its connection with quantum 
dots. In the equilibrium, Anderson model is exactly solv- 
able by the Bethe ansatz |2^, and it is a potential can- 
didate to be solved exactly out of equilibrium. It can be 
shown j2^ that, e. g., susceptibility in equilibrium can be 
expanded into the power series in U/T which converges 
absolutely for |J7/r| < oo and rapidly attains the asymp- 
totic form. It is natural to assume that such convergence 
properties remain true also out of equlibrium. Thus, if 
Eq. (^3|) holds, then it must be satisfied for every term 
in the {/-perturbation expansion of both sides, in anal- 
ogy with the standard linear response theory |28|. We 
checked Eq. (|2^) for the first few terms. In the first or- 
der in U the self-energy depends on the occupation of the 
levels at J7 = 0. The occupation (n) obtained from the 
analaytic continuation of the Matsubara Green function 
of the non-interacting system, 



Sf/)-C/C;°(r,r- 



(26) 
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does not in general agree with the exact non-equilibrium 
occupation of the same system. 



U 
2" 



J:f^^=-^UG-+{t,t) = ^-^lm 



(27) 



2 2'kT 



+ 1 



.£0 + M 



27rT 



27rT 



+ 1 



.eo - M 



2t:T 



where '^'{x) is the derivative of logarithm of gamma- 
function and we have taken fj, = fj.^ = —fiR > 0. We 
do not exclude, however, that for some special choice of 
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parameters Eq. (|2^) might hold. Indeed, from the above 
exphcit expressions in the first order we see that, except 
for the equihbrium point /i = 0, the two expressions agree 
also at the symmetric point eo = out of equilibrium, 
where = E^j'^j — U/2. Of course, one cannot rely 
merely on the first order result. It just indicates that Eq. 
( p3| ) could hold in the true symmetric point eo = —U/2 
where all the energy levels (/i, — /u, eo, cq + U) are symmet- 
ric with respect to zero and an additional particle-hole 
symmetry is present. For eo = — C//2 one treats the eod^d 
term as a perturbation, and the first order contribution 
to self-energy vanishes: S^^^ = S^^^^ — 0, thus the non- 
trivial check for the symmetric point can be provided by 
the second order in U. As can be shown, only one dia- 
gram. Fig. 1, gives a non- vanishing contribution up to 
this order. 




FIG. 1. First non-vanishing contribution to the self-energy 
of the symmetric Anderson model. 



(30) 



The first non-trivial correction to the c?-electron self- 
energy is given by the one-loop diagrams shown in Fig. 
2. 



(a) 



(b) 



FIG. 2. The leading contributions to the self-energy of 
d-electron interacting with phonons (dashed line). 

The diagram in Fig. 2(b) is proportional to G^i^iT, r+0) 
and thus reduces to the analysis of the first order correc- 
tion for the Anderson model (26). The diagram shown 
in Fig. 2(a), however, contains a non-trivial to and fx de- 
pendence and can be easily evaluated analytically. One 
has 

I]ff)(ic^„,Z/i) = -J^TY,n^p)Gddi^n+np), (31) 



Unfortunately, in the second order the calculation is 
significantly more cumbersome and we were not able to 
complete it and obtain an explicit analytical result. We 
hope, thus, that Eq, (|2^) may hold for the special choice 
of the value of eo in the Anderson model, leaving it as 
an open question for future studies. Note also that ac- 
cording to the Luttinger-Ward procedure for calculating 
the ground state energy of interacting Fermi gas j2^ the 
first order correction to the self-energy (which is due to 
instantanious self-interaction) cancels against the correc- 
tion to chemical potential. This result does not apply 
to our case, since the chemical potential is kept fixed in 
our problem, while the total number of particles, being 
in fact infinite, is allowed to vary. 

As another example, which allows to perform a reliable 
check of Eq. ( p3|) avoiding above difficulties, we consider 
model (H) with the dissipation on impurity. Namely, let 
Hint describes the effects of phonons (photons) emission 
and absorption on the impurity |3Q| , 



Hint = "fd^dcj), 



(28) 



where ^ is a real phonon field and 7 is the electron- 
phonon coupling constant. For simplicity we consider 
a toy model where the spectrum of phonons consists of a 
single mode, LiJ{){k) — loq > Q. The phonons are assumed 
to be in thermodynamic equilibrium initially, with the 
energy 



H 



ph 



(29) 



The field (f) is quantized in terms of the bosonic operator 
[b, 6+] = 1 as follows: 



where ftp is the bosonic Matsubara frequency and T>{t) 
— (T^(/)(t)0(O)) is the phonon Matsubara propagator. 



(32) 



The corresponding non-equilibrium retarded self-energy 
(Fig. 3) is given by 

/OO 7 

-D+-{u;,)G^+{u;^+u;)], (33) 

where the Keldysh correlators D are given in 

Appendix B. Substituting corresponding bare propaga- 
tors and evaluating (|3T|) , (|33| ) explicitly one can see that 
the only case when Eq. ( |23|) holds is the case of equilib- 
rium, /i = 0. 



FIG. 3. Leading contribution to the retarded self-energy 
of d-electron interacting with phonons, = E ^ -I- E , 
non-equilibrium case. 



IV. CONCLUSION 

To summurize, a generalization of the linear response 
theory to the non-equilibrium case was discussed. This 
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generalization is based on the certain analytic continua- 
tion in the additional complex parameter which plays a 
role of external bias. The complex chemical potential 
is formally introduced in the equilibrium density ma- 
trix and analytically continued to the real values after 
averaging over all the states is performed. Thermody- 
namic averages turn out to have non-analytic depen- 
dence on the complex bias, with the cut along the plane 
Im^i = Im/ifl, as a result of divergences in the sum- 
mation over infinite-particle states. It was shown that 
formula for the current is valid in the linear response 
limit for systems of the type (|^) and is valid also out of 
equilibrium if condition ( ^3| ) holds. This condition, which 
is trivially satisfied by the free systems, Hint = 0, is not 
satisfied in general by the interacting systems. This does 
not mean, however, that the very idea of obtaining non- 
equillibrium current by exploiting analytic properties of 
the thermodynamic averages in the space of complex cou- 
pling constants is hopeless. The fact that Eq. (|l]) holds 
for the free systems indicates that there are no general 
grounds for ruling out this idea. As the perturbative 
analysis indicates, the propagators obtained by an ana- 
lytic continuation from the Matsubara ones have similar 
structure and resemble real non-equilibrium propagators. 
Hence, it would be interesting to derive an analogue of 
Eq. based on the analytic continuation in the com- 
plex bias which would be true for the generic interact- 
ing systems as well. One may proceed, e. g., from the 
Lehmann spectral representation jSlj l to obtain a relation 
between the equilibrium partition function and the non- 
equilibrium current by controlling two complex parame- 
ters, frequency and chemical potential. Formula as it 
is can be used, nevertheless, to obtain the exact trans- 
port through Anderson impurity in the linear response 
limit Is^ , if the exact partition function is supplied by 
the Bethe ansatz technique. It would be interesting to 
understand also in the Green functions language why Eq. 
M) works for the quantum Hall bar with a constriction 
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APPENDIX A: FREE CASE 



1. Point-contact model 

To understand how the analytic continuation in the 
complex chemical potential works, we consider the fol- 
lowing simple point-contact model : 



H 



(Al) 



This model can be viewed also as a particular case of the 
quantum Hall bar with an impurity [Q at = 1, when the 
Luttinger liquid on edges reduces to free fermions, and 
thus it provides a direct check of Eq. (|l|) . By definition, 
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(A2) 



where pL = etp'^ipL is the density of charge. Opening the 
commutator above and taking the average, one obtains 



(A3) 



The latter expression involves equal time correlators at 
X — 0, which can be easily found by using Dyson equa- 
tions and the Keldysh technique ||T^. Eventually, one 
derives 



(A4) 



/i(i + r/4)2 



R 

{^^R. - ml) (A5) 



where fL,R — f{'^— ^J^L,R) are the Fermi functions, F = 
and is the retarded Green function of either of the 
disconnected leads (A = 0) at x = Pj, 



dp 



1 



27r cj + p + iO 



(A6) 



Partition function of (|AlJ) can be easily found. For the 
complex chemical potentials i^l and i/i^, according to 
(0) one has 



d 



log Zi„p(iAiL, iAifl) = i-vg^ 



5oL('^n)5ofl(Wn) 



— rCJoL(Wn)t/0_R(t^n) 
1 



2-kT F - 4 



^ F-f 4' 

„>max{ni,,fifi} 
„<min{ni,,Mjj} 



(A7) 



where 



1/T 



dre-"-(T.^i(«)(T)V'+^,)(0))o (A8) 



are the Matsubara Green functions for the disconnected 
leads (A = 0) at X = 0, e. g. 



dp 



1 



-sign(w„ - Ml), 



(A9) 
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and ujn are the fermionic Matsubara frequencies. When 



passing from the first line to the second in (A7) it was 
used that /il = 27rNT and fiR ^2nMT, where iV, M are 



integers. The infinite sum in (A?) needs to be properly 
regularized: 



E 



1 



1 



1 



1 >max[n£, 
T<min[fi£_,fiji] 



r + 4 



lim 



2mT r + 4 £^+0 e'^/^ + 1 



duj 



2ttT r + 4 ' 



(AlO) 



where contour C consists of two parts: one goes along 
Imw = max[/ii, /i/j] and an infinite semicircle above it, 
the other along Imcj = min[/ZL, /^i?] and a semicircle be- 
low it, both counter-clockwise. 

Considering /i^ — iin as one variable y, one has a func- 
tion \y\ defined on the imaginary y-axis. We have to 
continue it analytically to the real axis. Since it has a 
cut at ?/ = 0, one can define two separate analytic contin- 
uations, to the upper and lower half-planes: /+ — —iz, 
/_ = iz. Then, one continues from the positive imag- 
inary axis to the real axis, f+{x) = —ix, and from the 
negative imaginary axis, f-{x) ~ ix, and subtracts one 
from another to get f+ — f- — —2i{^L — M_r) fH]. Com- 
paring to (A5), one has 



j=.£87rzT^^log^. 



(All) 



Note that the prefactor (4 - r)/(4 T) in ( |AllD is not 
universal and depends on the regularization scheme. 



2. Resonant-level model 

Next, instead of the point-contact scattering we intro- 
duce a single dynamical impurity at a; = described by 
an additional state at the energy eg- Operator d"*" creates 
an electron on the impurity. The simplest model to write 
is a resonant-level model: 



(A12) 



m—L.R 



Note that the impurity can be effectively removed from 
the action by integrating out impurity degrees of free- 
dom, and one gets certain generalization of the point- 
contact model with time-depende nt coupling. Then, in 
the limit eo , A — > oo one recovers ( Al ) . A dynamical im- 
purity, however, suggests a solid way of regularizing the 
divergences encountered in Eq. (A7). 

Expression for the current was derived in Sec. |l|: 

\{Jl - Jr) = ^ / dLo[jL- fR]lmG^ 



J 



(A13) 



where for the retarded Green function one has in this 
case an exact expression 



G: 



1 



dd 



uj — eo + iV 



(A14) 



Substituting (A14) into (A13) one obtains 



-T Im 



27rr 



+ 1 



2ttT * 2ttT 



2-kT 



(A15) 



where ^' is the standard psi (digamma) function (com- 
pare to Eq. dH)). 

The partition function involves gaussian Feynmann in- 
tegral and can be obtained immediately: 



d 



E 



iujn - eo - r(t/oL + Gor) 
i 



iuJn + i max[^L, Aifi] - eo + «r 



iuJn - i inm[fiL,^ii] + eo + iV 



(A16) 



where, while passing from the first line to the second 
above, Eq. (A9) was used and the fact that ifiL.R have 



the form of bosonic Matsubara frequencies, which allows 
to bring ijJL^R into the denominators. Expression ( Al(^ ), 
viewed as a function of two complex variables /il and 
fiR, has a cut at /ii = /^_r. The analytic continuation 
is performed from the 2D plane of purely imaginary fi^ 
and fiR with the cut a,t = ^r. Once one continues 
analytically from the domain fi^ > I^r, and once from 
Mi? > Mi to get two functions. Subtracting one from 
anoth er an d substituting i/i — > /u, one gets by comparing 
with ( |A15| ) 



J 



3. Kondo model 



(A17) 



Consider an impurity spin S = which couples L and 
R leads by the s — d exchange (Kondo) interaction pSl : 



Hk^Ho 



3 

E E 

A=l ab=L,R 



Jfib+M'TL'M0)S\ (A18) 



where are Pauli matrices and J^^ are coupling con- 
stants. Electrons here have additional spin index a. As 
it was shown in following the ideas of , under the 
assumptions that ^ = J^-^ = 0, Jx = Jy = J±, Jf ^ = 



7 



jRR^ jLR ^ jRL^ jLL ^ jRR ^ ^TX model ( |AT8| ) is 

equivalent to a quadratic model with the Hamiltonian 
ilK= J2 ^fdx^P+d^^Pm + JiS{x){4>++^l;i){d+-d) 

(A19) 



m=l,2 

+ J2S{x){^+ -^P2)id+ +d), 



where Ji ~ Jf-^, J2 ^ Jf-'^, and {d+,c?} = 0. Indexes 1 
and 2 correspond to the spin-flavor and flavor channels of 
the original model (AI8) [p6|, while only second channel 



has now non-zero chemical potential /i if one sets initially 
fj-L = —fJ-R = §• When calculating Zimp it is convenient 
to change variables in the Fcnmann integral from d, to 
the real and imaginary parts a — {d + d~^)/ -s/S, b — i{d — 
(i+)/-\/2, where a and b are classical Majorana fermions. 
After integrating out ?/'i,2 the action factorizes as = 
S{a, j|, fi) + S{b, Jf ), and only the /i-dependent part is 
of interest to us. Since S is gaussian, Zimv can be easily 
computed. The result coincides with Eq. (Al£) at eg = 
0, — ~M_R- Thus, Kondo interactions reduce to the 
case of resonant-level model studied above for the special 
point in the couplings space (an analogue of Toulouse 
point for the single-channel Kondo problem). 



APPENDIX B: GREEN FUNCTIONS 

We collect here the definitions of different Green func- 
tions used throughout this Letter, following the notations 
of Ref. All the Green functions are taken at a; = 0. 
The subscript R for "right" should not be confused with 
superscript R for "retarded". First, consider the point- 
contact model. Define: 



i{^tmR{t))^Gl+iO) 



(Bl) 
(B2) 



Dyson equations yield: 



GZi = >^9R-Glt ^9r^GU (B3a) 



Grx — G_Rfl ^ ^9l^G^'^ 



(B3b) 



where 



dte-*(fiV'L(i?)(t)^+(^)(0)) (B7) 



(B8) 



are Green functions of free (decoupled, A = 0) leads, Tt 
stands for anti time ordering operation, and by capital 
letters Gll,Grr are denoted corresponding full propa- 
gators of the interacting system (A 7^ 0). The latter can 
be found from the system of four coupled Dyson equa- 
tions. The solution reads, e. g.. 



+ _ 9Ll±>Md9Rl_ 



|l-A2gf5gl2 



r, + + J- \2„_R„A^-- 

++ _ 9l + ^ 9l 9l 9r 



|l-A2gf5g'^ 



(B9) 
(BIO) 



For the resonant-level model there are very similar defi- 
nitions. For example. 



t{d+{t)Mt))^G-+{0) 



z(V'+(t)d(0)=G£+(0) 



(Bll) 
(B12) 



Dyson equations express these correlators through the 
full propagator of the impurity: 



GdL - ^9l ^dd ^ ^9l G^^ 
GlI = -^^dd 9l^ - ^G^i^gj^^ 



(B13a) 
(B13b) 



Four Green functions 4 
lated to each other by 



, H — , — -f in general are re- 
(B14) 



G++ + G— = G+- +G-+ 
and to the retarded (advanced) Green functions 

G^ ^G— -G-+, (B15) 



G"^ = G-- - G^ 



(B16) 



The impurity Green functions of an isolated non- 
interacting impurity read: 



G 



-+{o) 

dd 



iiTd{u! — Co) 



(B17) 



^9LiR)('') = / rfte*"*(^^(;j)(0)^L(ii)(0> 
J —00 

= f{uj-tiLiR)), (B4) 



(B5) 



*5Z(fl,)(^) = / dte^'^*(TtVL(i?)(i)V'2(^)(0)) (B6) 



'9t[n)(^) = / dte*'^*(^L(i^)(0^^(;^.)(o)) 

•/ — 00 

= [l-/(^-AiL(i?.))] 



G 



+-(0) 



dd 



G 



fl(0) _ 
dd 



-ITtS^LU — Eq) 



1 

LO — Eq + iO 



(BIS) 



(B19) 



For the resonant-level model without on-site interactions 
one can write four Dyson equations for the impurity cor- 
relators and easily solve them. Relevant for us is the 
retarded Green function. 



r<R — 

^dd — 



G 



RiO) 
dd 



gT^- 



(B20) 
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where T,^ = T, +E + is the retarded self-energy, given 
in the free case by 



(B21) 



For the calculation of interacting Green functions one 
needs also the following zeroth-order propagators: 



Gdd 



Gdd - 



_ '-'dd 




9r^) 


|1 


-^'GTi9S + 9^W 




^dd 


^>^'GTGTi9t- + 




|1 






r — (0) 
_ ^dd 


,X^GTGTi9t^ + 


9lt) 


|1 







(B22) 



(B23) 



(B24) 



Zeroth order correlators of the phonon field are defined 
as 



(B25) 
(B26) 



Explicitly, 



D-+{uj) = -iTTUJo[N6{uj - ujo) + (1 + N)6{uj + ujq)] (B27) 
D+-{u;) = -i7:LUo[N6{LJ + wq) + (1 + N)S{u; - luq)] (B28) 
1 1 



with 



2 



uj — ujQ + iO uj + luq + iO 



N = 



1 



(B29) 



(B30) 
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